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Superconducting states of an anisortopic s-wave superconductor on a Mobius strip are studied 
numerically based on the Ginzburg-bandau theory and the Bogoliubov-de Gennes theory. In both, 
the equations are solved numerically on discitized lattice and the nonlinearity and the self-consistency 
are fully taken into account. First, we study the superconducting states on the Mobius strip in the 
presence of the Aharonov-Bohm flux threading the ring by employing the Ginzburg-Landau theory, 
and confirm the phase diagram previously proposed by Hayashi and Ebisawa [J. Phys. Soc. Jpn. 
70, 3495 (2002)]. The metastable states as well as the equilibrium energy state are studied and 
the nonequiriblium processes when the magnetic field is varied at a fixed temperature are discussed. 
Next, we study the microscopic superconducting states on the Mobius strip based on the Bogoliubov- 
de Gennes theory, especially focusing on the state with a real-space node in the superconducting 
gap, which is expected to appear when the flux threading the ring is half the superconducting flux 
quantum. The local density of states in this nodal state is calculated in detail and the existence of 
the zero-energy bound states is shown. 

PACS numbers: 74.20.De, 74.78.Na 



I. INTRODUCTION 

The realization of crystals with unusual shapes, e.g., 
ring, cylinder etc., by Tanda et al^2£ has stimulated 
renewed interest in the effects of the system geometry 
on the physical properties. Especially, the synthesis 
of Mobius strip made of transition metal calcogenides 
(NbSe3, TaS3 etc.)^ opens new possibility to examine the 
physical properties of superconductivity or charge den- 
sity wave in topologically nontrivial spaces. 

Recently several groups have studied physical systems 
on Mobius strips. Hayashi and Ebisawa 5 studied s-wave 
superconducting (SC) states on a Mobius strip based 
on the Ginzburg-Landau (GL) theory and found that 
the Little-Parks oscillation, which is characteristic to the 
ring-shaped superconductor, is modified for the Mobius 
strip and a new state, which does not appear for ordi- 
nary ring, shows up when the number of the magnetic 
flux quanta threading the ring is close to a half odd in- 
teger. Yakubo, Avishai and Cohen& have studied the 
spectral properties of the metallic Mobius strip with im- 
purities and clarified statistical characteristics of the fluc- 
tuation of the persistent current as a function of the mag- 
netix flux threding the ring. The persistent current in a 
more simplified version of the Mobius strip has also been 
studied by Mila, Stafford and Caponfi. Wakabayashi 
and Harigaya& have studeid the Mobius strip made of 
a nanographite ribbon, and the effects of Mobius geome- 
try on the edge localized states, which is peculiar to the 
graphite ribbon, has been clarified. A study from a more 
fundamental point of view can be found in the paper by 



Kaneda and Okabe^ where the Ising model on Mobius 
strip and its domain wall structures are studied. 

The main result of Ref. is that if the magnetic flux 
threading the ring is close to a half odd integer times 
the flux quantum 4>o = ftc/(2e) (h, c, e being the Planck 
constant, the speed of right and the electron charge, re- 
spectively), a novel SC state appears. This state has a 
real-space node in SC gap along the strip: namely, the 
gap tends to zero along the line located in the middle of 
the strip. Throughout this paper we call this state the 
11 nodal state 11 . It has been shown that the free energy 
of this state can be lower than that of uniformly gapped 
state, which is known to be the most stable state in case 
of the ordinary SC rings. 

This paper extends the previous study in the following 
two points: 

a) In Ref. |f| we have shown that the free energy of the 
nodal state is lower than other likely states. How- 
ever, there is no evidence that it is the most stable. 
Although it is difficult to examine all possible local- 
minimum states of the nonlinear GL free energy, in 
this paper we try to give a more convincing evi- 
dence by resorting to a numerical method. We per- 
form numerical minimization of the GL free energy 
and find as many local-minimum states as possible 
and reexamined the phase diagram of Ref. H A 
similar method is previously employed, for exam- 
ple, in Refs. ^ andlTTl in studying SC disks etc. 

b) The electronic states, which are not treated in Ref. 
0, are studied in terms of the Bogoliubov-de Gennes 
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FIG. 2: Phase diagram of a sperconducting Mobius strip 
FIG. 1: (a) Geometry of SC Mobius strip where an Aharonov- in the presence of a AB flux with (a) < r± and (b) 

Bohm flux is threading the ring (b) The developed figure of HUL > r± 
(a). 2v ^ 



(BdG) theory. We numerically solve the BdG equa- 
tion on a lattice selfconsistently, where electrons are 
treated by the tight-binding approximation. Based 
on the solution we study the local density of states 
in the nodal state, which may be observed, for ex- 
ample, by scanning tunneling microscope measure- 
ment. 

This paper is organized as follows: In Sec. ^ we re- 
view the behavior of a SC Mobius strip in a magnetic 
field, presented in Ref. la. In Sec. IIIII the model, method 
and results of the numerical analyses of GL theory are 
presented. In Sec. IIVI studies based on the BdG theory 
for the electronic properties of the nodal state are pre- 
sented. In Sec. E] discussions on the results and their 
relation to experiments are given. Sec. IVII is devoted to 
summary. 

II. A SUPERCONDUCTING MOBIUS STRIP IN 
A MAGNETIC FIELD 

Here we summarize the results obtained in Ref. 
In that paper, the behavior of a SC Mobius strip in an 
external magnetic field is studied based on the GL theory. 

We consider the Mobius strip made of a superconduc- 
tor as shown in Fig. 2] (a). The magnetic field is as- 
sumed to be threading the ring in a form of the so-called 
Aharonov-Bohm (AB) flux as indicated in the figure by 
a bold arrow, which gives rise to a non-zero vector po- 
tential on the strip although the magnetic field is vanish- 



ing on the strip. We further assume that the coherence 
length along the strip and that perpendicular to the strip, 
which are respectively denoted as £m and are different. 
Throughout this paper we set x-axis and y-axis as indi- 
cated in the figure. The developed figure of the Mobius 
strip is given in Fig. H(b). We denote the circumference 
and the width by L and W, respectively. 

In this paper the effects of the bending, or the non- 
zero curvature, caused in the strip by the Mobius geom- 
etry, are neglected. Although these may be important 
in the case of unconventional pairing, such as p-wave 
superconductors 12 , they may be neglected for s-wave su- 
perconductors, treated in this paper. 

The equiriblium states of SC Mobius strip depends on 
the magnetic flux $ threading the ring. The SC Mobius 
strip shows the so-called Little-Parks oscillation of the 
transition temperature with a period <po, as we naively 
expect from the analogy to the ordinary rings. However, 
it turned out that the oscillation can be appreciably mod- 
ified in SC Mobius strip depending on the strength of the 
anisotropy of the coherence lengths. 

The analysis of Ref. 5 : shows that there are two impor- 
tant parameters, which are given by 



where £_i_(0) and £y(0) are the coherence lengths at ab- 
solute zero temperature(T = 0). 

It has been shown that when the condition — 2=7-11 < r± 

2V3 1 

is satisfied, the phase diagram in a magnetic field behaves 
like the one shown in Fig. [21 (a), which is basically the 
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FIG. 3: (a) Real part, (b) imaginray part and (c) amplitude 
of the order parameter in the nodal state. 



same as in the case of ordinary ring. Here T c is the tran- 
sition temperature in the bulk and the index n indicated 
in the figure denotes the the winding number of the phase 
as we go around the ring once along a trajectory parallel 
to the edge. When 



2 v / 3 



(2) 



is satisfied, a state characteristic to the Mobius geometry 
appears when the number of the flux is close to a half 
odd integer, as shown in Fig. [2] (b) by hatched regions. 
These states are indexed with a half odd n and the phase 
of the order parameter changes by an odd number times 
7r as we go around ring once. The spatial dependence of 
the order parameter is shown in Fig. |3J where the real 
part, the imaginary part and the amplitude are shown 
for n = 1/2 case. It is clear that the order parameter has 
a real-space node in the middle of the strip. Therefore 



we call these states the "nodal states" . The nodal states 
can exist above t\ shown in Fig. 121(b). Here t\ is given 
by 



ti = l- 



3tt 



2 r\ 



(3) 



as one can see from the results in Ref. 0- 

The results of Ref. |j| are obtained by comparing the 
free energies of the several possible states, which are cho- 
sen empirically. Therefore it is not easy to say that there 
is no states with lower free energies. Of course, it is im- 
possible to investigate all possible order parameter con- 
figurations. However more reliable analysis, which can 
cover wider range of the configuration space is required. 
In this paper, to fulfill this requirement, we perform nu- 
merical study on GL theory, which is given in Sec. EH 

The analysis of Ref. is limited to the phenomeno- 
logical one. Since the nodal state is also interesting from 
electronic points of view, more microscopic study is re- 
quired. We solve BdG equations on the Mobius strip 
numerically and clarify the electronic bound state near 
the node of the nodal state. This will be given in Sec. 

EH 



III. NUMERICL STUDY BASED ON 
GINZBURG-LANDAU THOERY 

In this section, we study the SC state on the Mobius 
strip numerically using the GL theory. Here we employ 
non-linear optimization method (quasi-Newton method) 
to find the local minimum state of the GL free energy^. 

GL free energy F of our system is given as follows: 



F = di 



d z r 



2m* 



8 .2tt, . 

75--*-7-4e V 

OX (j)o 



ly-% Ay]l 



+ a (t-l)|Vf + f|V>| 4 



(4) 



Here m* , <fio = hc/(2e) are the mass of a Cooper pair 
(twice the electron mass) and the magnetic flux quan- 
tum, respectively. The thickness of the strip d± is much 
smaller than the superconducting coherence length and 
the strip can be treated as two-dimensional. ip(f) and 
A(f) are the SC order parameter and the vector poten- 
tial, respectively, t = T/T c is the reduced temperature. 
The AB flux is incorporated by taking A — $>/Le x , where 
e x is the unit vector in x-direction. ao and (3 are positive 
constants, and 7 = is the anisotropy parameter. 

For numerical calculations we introduce the lattice ver- 
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sion of F as 



F =F< 



L j=i fc=i 



c.c. 



AT X JV„-1 

+t 2 E £ f |^0",*)-^0',*+i) 

3=1 fc=l 

+EEi( < - 1 )i^^)i 2 + 5i^*)i 4 



(5) 



»=i i=i <■ 

where o x = 2TrdA x /4>o, £ is given by £n(0)/d, where 
£|| (0) = y / h 2 /2m*a a and d is the lattice spacing. Here 
the lattice is assumed to be a square one. Note that 
N x d = L and N y d = W. F is V Q a 2 /P where V = 
d 2 x d±. The order parameter ip is normalized so that 
ijj — > 1 as t — > 0. 

We use the open boundary condition in y-direction and 
Mobius boundary condition in .T-direction: namely, we 
put 



^(JV a + l,fc) = ^(l,JV w + l-*), 



(6) 



(1 < k < N y ) in the first two lines of Eq. ©. 

The order parameter ip(j, k) minimizing F is obtained 
by solving the equation, 



dF 



dip* (j, k) 



0. 



(7) 



which yields, 

f 2 {-'0C? + 1, k)e-^ - $(j - 1, fc)e i5 « 

- 7 2 ^(j, k + 1) - 7 2 ^(J, fc - 1) + 2(1 + 7 2 )^(i, fc) 

+ (t- l)^(j,fc) + |^(j,fc)| 2 VXj^) = 0. (8) 

We obtained the solution of the Eq. (JSJ in terms of the 
nonlinear optimization of the free energy Eq. (JSJ) . In this 
paper we especially utilized the so-called quasi-Newton 
method 13 . Like other methods of nonlinear optimization, 
quasi-Newton method starts from an initial value of ip 
and changes it so that the free energy becomes lower until 
finally we reach the local minimum. Therefore which 
local minimum we reach depends on the initial value of 
ip. In this paper we randomly chose the initial values 
and performed optimization as many times as possible. 
Then we obtained several local-minimum states. These 
procedure is the same as the one adopted in Refs. ^] and 

a 

The results are as follows. The system size we used is 
N x = 10, N y = 10, Co = l-5d and £ ± = 1.2d. 
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FIG. 4: The free energy of the most stable state as a function 
of T and <!>. Around $ ~ (po/2 a structure corresponding to 
the nodal state can be seen. The vertical axis is the logarithm 
of the \F\ and the region with \F\ < 10~ 5 near T c has been 
cut off. The free energy is measured in the unit of Fq. 
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FIG. 5: The free energy of the metastable states as a function 
of $ for (a) t = 0.78, (b) t = 0.5, and (c) t = 0.1. The 
branchs corresponding to the nodal states are indicated by 
bold arrows in every graph. It becomes the most stable state 
when t = 0.78 and 0.5 for a certain region of "I". When t = 0.1 
there is no stable region for the nodal state. 
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First we discuss the phase diagram. The numerically 
obtained free energy is shown in Fig. fusing the log scale: 
we have shown only the region with F < and log \ F\ 
is plotted as a function of $ and T. We can see from 
this figure that the structure corresponding to the nodal 
state appears near <I> = </>o/2 (shown by bold arrows). 
Actually, the order parameter in this region behaves like 
that shown in Fig. [21 

In Fig. we have depicted the free energies of the 
metastable states at three different temperatures ((a) t = 
0.78, (b) t = 0.5 and (c) t = 0.1). We clearly see three 
distinct series of states in each graph. The branch which 
starts from (f> = corresponds to the state with a uniform 
order parameter A = (const.) and that ends at 4> = 00 
corresponds to the state with A = (const.) x e l27TX / L ^ 
Between these two branches, a branch corresponding to 
the nodal state can be seen. At t — 0.78 and t = 0.5, 
the nodal state has the lowest free energy near 4> ~ 4>o/2. 
However, at t — 0.1, although the nodal state is still 
a metastable state, there is no region where the nodal 
state is the most stable. These features agrees with the 
prediction of Ref. 0. 



IV. NUMERICL STUDY BASED ON 
BOGOLIUBOV-DE GENNES THOERY 

In order to study the s-wave SC state on a Mobius 
strip microscopically, we treat a tight-binding model on 
a square lattice with attractive on-site interactions. The 
Hamiltonian of the system is given as 



with A 



H 



-t S^fp^A r 
3" 

~ty ^^( C j+y,cr C 3,^ C j,a c 3+y,v 
3° 



p -i0x„t _ ) 
e L -j,a L '3+ x , a I 



(9) 



where c JCT is the annihilation operator of electron at a site 
j with spin cr(=T, I), and V(> 0) and [i are the attrac- 
tive interaction and the chemical potential, respectively. 
Here j = (j x ,j y ) (1 < j x < N x ,l < j y < N y ) num- 
bers the sites, where N x and N y are the numbers of sites 
along the x- and y-direction, respectively, and x = (1, 0), 
y = (0,1). n,j a = Cj^Cja is the electron number oper- 
ator. The transfer integrals for x- and y-directions are 
denoted as t x and t Vl respectively, and the Peierls phase 
4>x = (tv/N x )(^/4>o) represents the effect of the AB flux 
$ threading the Mobius strip. To realize the Mobius ge- 
ometry we impose, for x-direction, the condition, 

C (N x + l,j y ) : a = C{l,Ny-j y +l),<Ti (10) 

while for the y-direction we assume the open boundary 
condition. 

The interaction term is decoupled within a mean-field 
approximation as 



n n n jL -> A j c] i c] T + A*c n c jL - |A 3 f 



j = \Cj]Cji) being the SC order parameter. Then 

(12) 



the mean-field Hamiltonian is written as 



where 



with 



hjk = 

W jk = 
F jk = 



W 3k 



Fjk 

-w jk 



it 

C 3l 



By solving the following BdG equation 



I 



Uln 
Vln 



— E n 



(13) 



(14) 



(15) 



we can obtain the energy eigenvalues E n and the corre- 



sponding eigenfunctions (v,j 



, ) , where n is numbering 



the states. The unitary transformation using (uj n ,Vj n ) 
diagonalizes Hmfa, and the SC order parameter Aj can 
be written in terms of E n and (uj n ,Vj n ). These consti- 
tute the self-consistency equations which will be solved 
numerically. In the following we take the parameters as, 



t x = 1.0, t y = 0.49 

fi = 0, T = 0.22, y = 0.25. 



(16) 



First we estimate the correlation lengths, £n (0) and 
C-l(0), by calculating the anomalous Green's function 
F(J) = {ck[Ck+j]) at T — without applying the mag- 
netic flux, and the results are depicted in Fig. (This 
calculation has been carried out in a larger system with 
100 x 100 sites using periodic boundary condition for both 
x- and y-directions.) From the results in Fig. [S]and the 
relation F(j) oc exp(— |j|/£), £||(0) an( i are esti- 

mated as 



MO) = 2.37d, ^ ± (0) = 0.64d, 



(17) 



(11) 



where d is the lattice spacing. 

The calculation in the Mobius geometry is carried out 
for a system with N x — 13 and N y = 14. Putting 
L = 13 d and W = 14 we obtain r|| = 0.18 and 
r± = 0.046. These values satisfies the condition, Eq. 
Furthermore, t\ is estimated to be —13, which means that 
the nodal state is stable at all temperature range down 
to OK. In passing, the bulk T c of this system is 0.358 (es- 
timated numerically in the system with 100 x 100 sites), 
which coincides with that in 13 x 14 system with $ = 
within a numerical accuracy. 

In this calculation, we limited ourselves to the case 
of the nodal state at <f> = (f>o/2. Since this microscopic 
calculation needs more time to obtain good convergence 
as compared to GL calculation, we have selected for the 
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FIG. 6: The anomalous correlation function (ckiCk+jt) f° r 
the system with the parameters in Eq. 11160 . (a) and (b) cor- 
respond to the correlation in the x- and y-direction, respec- 
tively. The equations in the figures are the approximating 
lines to the exponential decay. 




FIG. 7: The local density of state in the nodal state on the 
Mobius strip 




initial order parameter the solution obtained by the GL 
analysis, such as the one depicted in Fig. [3] After the 
calculation, the behavior of the order parameter has not 
changed so much and we consider that the iteration con- 
verged to the nodal state. The local density of states 
(LDOS) is calculated from the equation, 

N(j,E) = --ImY . (18) 

n 

where T is the broadening of the single energy level, in- 
troduced to simulate the actual experiment. The result 
is shown in Fig. [3 with T = 0.03. Each line of Fig. 
corresponds to the LDOS at j — (l,j y ) where j y num- 
bers the chain from the edge of the strip. (The LDOS is 
independent of j x .) Because of the inversion symmetry 
with respect to the center, we depicted LDOS only for 
1 < jy < 7, where j y — 1 and j y — 7 correspond to the 
outermost and the innermost chain, respectively. 

In Fig. a well-developed gap behavior can be seen 
for 1 < j y < 5 and the bound states are formed in the 
chains j y = 6 and 7. These bound states originates from 
the node. Since the node can be regarded as the sponta- 
neously formed 7r-junction (as one can see from Fig. |3 
the phase of the order parameter changes by it in crossing 
the node), the bound state energy must be zero. In Fig. 
{7\ the deviation from zero energy is seen, which may be 
due to the finite size effect. This point will be discussed 
in detail in a separate paper—. 
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FIG. 8: The transition of the states in the field sweep exper- 
iment for (a) t = 0.5 and (b) t — 0.78. The nodal state is the 
most stable in the shaded region. See text for the details. 



V. DISCUSSION 

In this paper, by using the numerical methods, we have 
confirmed that the predictions of Ref. are satisfied 
within both GL and BdG level. Until now, we have lim- 
ited our discussion to purely two-dimensional cases. This 
is not the case for actually synthesized Mobius crystals^. 
Here we discuss the effects of the finite thickness of the 
strip. If the Mobius strip is thicker than the SC coher- 
ence length perpendicular to the strip surface, the nodal 
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state may no longer be stable and the node, which is 
obtained under the assumption that the gap is uniform 
in the direction of the thickness, may become a vortex 
line, embedded inside the strip. Then there is no gapless 
region on the surface. In this case, we can not observe 
the node by only measuring the surface density of states 
by scanning tunneling microscope. The observation of 
the nodal state may be possible only by measuring the 
variation of the gap or magnetization as a function of 
the magnetic flux. More precise analysis on the thicker 
Mobius strips is left for the future study. 

In Sec. IIIII we have calculated the free energies not 
only of the equilibrium state but also of the metastable 
states. The metastable states do not appear in the 
thermodynamic equilibrium, though they play important 
roles in actual experiments in which one sweeps the mag- 
netic flvrsii^iSiiLiS. As shown by the previous studies, 
when we change the magnetic flux the state of the sys- 
tem remains in one of the metastable states and it does 
not switches to another lower-energy state until the ini- 
tial state finally becomes unstable. Although there may 
be effect of thermally assisted tunneling, such effects are 
limited to the very vicinity of the critical temperature. 
From these we point out a possibility that the nodal state 
may not appear at all in the field sweep experiment even 
though it can be a true equilibrium state for some values 
of the flux. This is understood from the Fig. [S] (a) , where 
the free energy for t = 0.5 and the time evolution of the 
system under a field sweep process are shown. The bold 
curves show the free energies of the metastable states and 
the equilibrium state at t = 0.5. The branch correspond- 
ing to the nodal state is indicated by a bold arrow. The 
dashed and dotted curves show the time evolution of the 



system in up- and down-sweep experiment, respectively, 
where it is assumed that the transition between differ- 
ent branches are prohibited by the energy barrier except 
when the system comes to the end of a branch. In this 
case, although the nodal state can be a true equilibrium 
state, it does not appear during the field sweep process. 
In contrast to thie, the nodal state appears during the 
field sweep process at t = 0.78, as shown in Fig. (b) 
(Only up-sweep process is indicated). Therefore the ob- 
servability of the nodal state in the field-sweep experi- 
ment may be further limited to the region in the vicinity 
of the critical temperature. More precise numerical sim- 
ulation is required to clarify these dynamical processes of 
this system quantitatively, which is left for future studies. 



VI. SUMMARY 

In this paper, we have investigated the superconduct- 
ing states on a Mobius strip in terms of the numerical 
analyses based on Ginzburg-Landau and Bogoliubov-de 
Gennes theory. It has been shown that in the Mobius ge- 
ometry a novel nodal state can appear both in equilibrium 
and metastable states. The experimental observability of 
the nodal states is discussed based on the findings. 
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